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^ As in such 2D systems as the XY model, topological (or singularity point) defects are thought 

to play a crucial role in the phase transitions of 3D spin systems. In double exchange (DE) 
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ferromagnets, the conduction electrons are strongly coupled with core spins through Hund's rule, 
^ and in the presence of non-trivial spin texture, acquire a Berry phase contribution to the anomalous 

^ Hall effect (AHE). We combine Hall and magnetization data on CrC>2 with a thermodynamical 

o 

{Nl scaling hypothesis to confirm that the critical behavior of the topological-spin-defect density is 
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consistent with that of the heat capacity This analysis is the first experimental confirmation of 
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£h the topological character of critical fluctuations. 
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Topologically stable defects in materials often play a key role in phase transition phenom- 
ena. Topological defects (or vortices) are known to drive phase transitions in a large class 
of 2D systems, including superfluid films, 2D crystals, and planar-spin (XY) models [1, 2]. 
In 2D systems, singular excitations (topological defects) can be separated from nonsingular 
ones (spin waves) and treated analytically, leading to the conclusion that topological defects 
are both necessary and sufficient to describe the critical properties. Because the role of 
topological defects in 3D systems has not been investigated analytically, we must rely on 
Monte Carlo (MC) calculations [3-5] for evidence that topological defects play a crucial role 
in the phase transitions of, e.g. 3D Heisenberg ferromagnets. In 1989, Lau and Dasgupta 
(LD)[4] performed MC calculations on the classical 3D Heisenberg model, both allowing 
and suppressing topological defects. Without such defects, the system remains ordered at 
all temperatures, indicating that topological excitations play a central role in the phase tran- 
sition. Allowing topological defects, they concluded that d(n)/dT ~ \t\~^, where (n) is the 
defect density, t = (T — Tc)/Tc, ip ~ 0.65 and suggested that (n) is a disorder parameter. 
Several years later, however, Holm and Janke (HJ)[3] repeated the MC simulation using 
larger lattices. They reasoned that quasi-local defects should have energy-like character and, 
contrary to LD's conclusion, found weak evidence that d{n)/dT ~ \t\~ a ] i.e., controlled by 
the heat-capacity exponent a. While the MC evidence is strong, it is very important to 
establish the existence of such defects experimentally and to study how they are related to 
the ferromagnetic-paramagnetic phase transition. 

The topological defects detected by LD produce a non-trivial spin texture that, in the 
presence of strong coupling between the conduction-electron and core spins, contributes to 
the Hall effect via the induced Berry phase [6-10] . As this contribution requires the presence 
of non-zero magnetization (and spin-orbit coupling) we denote it here as "the anomalous 
Hall effect (AHE)," even though it is not simply proportional to the magnetization. We 
assert that in double-exchange (DE) systems the strong coupling between the spin of the 
atomic core and that of the mobile electron make the AHE a sensitive probe of the nature 
of topological spin defects in 3D systems. 

According to recent insights into the band structure of Cr02 by Korotin et al.[ll], Cr02 
can be regarded as a self-doped DE system due to strong hybridization of d yz and d zx orbitals 
with oxygen p orbitals and Hund's-rule coupling of these itinerant (i-electrons with the 
localized d TV core. With one <i-electron localized and one itinerant Cr 4+ is similar to Mn ions 
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in some doped manganites that show colossal magnetoresistance (CMR). Unlike these CMR 
manganites, &O2 has no strong Jahn- Teller effects, no disorder and no formally different 
valence [11]. Except for the origin of ferromagnetism, CrC>2 is more like a conventional 
ferromagnetic metal such as Ni and therefore can be an appropriate material in which to 
study pure DE ferromagnet properties. 

CrC>2 films were grown by ambient pressure chemical vapor deposition on rutile (bet 
structure) Ti02 (100) substrates following the recipe developed by Ishibashi et al.[12] The 
saturation moment of this material at low temperature is close to 2 jUg/f.u., suggesting 
that the sample is stoichiometric. The sample reported here is ~ 2300 A thick and was 
patterned by photolithography and Ar-ion beam milling. In order to measure resistance 
and Hall voltages along the inequivalent in-plane directions, the sample was patterned as 
shown in the lower inset of Fig. 1; that is, the current path is "L" shaped, the bars aligned 
along 6((010))and c— axes ((001)) respectively. The width of the current path is 100 /itn 
and the distance between two voltage leads for the resistance measurement is 500 pm. The 
structural quality of the film was checked by x-ray diffraction; we could find no chromium 
oxide signals besides the correct phase of CrC>2(200). The off-specular reflections show 2- fold 
symmetry confirming epitaxial growth. The residual resistivity ratio (RRR = P^qqk/ Pwk) * s 
another criterion of quality for a metallic sample and is about 66 for c— axis and 20 for b— axis 
(Fig. 1), similar to earlier results [13]. The resistivity in both directions shows an inflection 
point around T = 380 K corresponding to Tc- As shown in the inset of Fig. 1, using the 
Fisher-Langer model for dp/dT [14], we determined Tc = 383.8 K and critical exponent 
a= -0.17; the latter agrees with the value from the scaling relation a + 2(3 + 7 = 2, with 
j3 = 0.37, 7 = 1.43 for CrC>2[15]. This agreement of the Fisher-Langer exponent with critical 
exponents from the magnetization measurements confirm that CrC>2 is a 3D Heisenberg- 
like ferromagnet. We note that a previous scaling analysis of a poly crystalline sample gave 
totally different exponent (a = —0.37) [16]. It appears that the Fisher-Langer relation is 
valid for this class of ferromagnet s despite the unconventional ferromagnetic mechanism and 
relatively low (but metallic) conductivity. 

In Fig. 2, we plot the field dependent Hall resistivity at various temperatures from 10 to 
480 K. The Hall resistivity of magnetic materials is usually expressed by 

p xy = E xy /j x = R B + R s p M, (1) 
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with R the ordinary and Rs the anomalous Hall coefficients. Throughout our experiment, 
magnetic fields were applied perpendicular to the film plane, (if//(100)) and therefore the 
demagnetization factor can be regarded as unity. As one can see, both directions of the Hall 
resistivity (open circles and lines) are identical within an experimental error, in agreement 
with Onsager's principle that p = —p yx regardless of crystal orientation. As noted above, 
we term p xy = p xy — R B the anomalous Hall resistivity, or AHE. Generally, the low tem- 
perature Hall effect arises only from the ordinary (OHE) if there is no intrinsic spin disorder 
or frustration [9] . The OHE is hole-like at low temperatures, in agreement with earlier work 
[13] but appears to have a kink in the low temperature p xy curves at around B = 1 T, prob- 
ably not due to AHE but rather OHE from the complicated multiband structure[17]. Above 
100 K, the AHE gradually arises and the p xy curves become large and negative, indicative 
of an electron-like AHE. With increasing temperature, the magnitude of p xy maximizes just 
below Tc(=384 K) with our achievable field, and gradually decreases. This behavior is qual- 
itatively quite similar to that of CMR systems and should be explained by the Berry phase 
theory [6, 10] rather than conventional spin scattering mechanisms such as side jump and 
skew scattering[18]. 

According to the MC simulations [3-6], topological defects are dilute at low temperatures 
with positive Q = +1 (hedgehogs) and negative Q = — l(anti- hedgehogs) separated by the 
lattice spacing a and strongly bound into pairs (or strings). Each pair acts as a lattice- 
spacing-size magnetic dipole. The temperature dependence of the topological spin-defect 
pair density n at low temperatures is given by (n) oc exp(— Ec/hsT), where Eq is the 
core energy of a pair. Well below Tq, the defect pairs are dilute and behave independently. 
When spin-orbit coupling X so is taken into account, a gauge field (b z ) appears due to a 
slight imbalance between the numbers of pair dipoles oriented parallel and antiparallel to 
magnetization M: 



Here, $ = hc/e = 4.1357 x 10~ 15 Wb, and n e /(~ 0.5) is the number of carriers per unit cell. 
We neglect higher topological charge defects (\Q\ > 2). Because X so is small compared to 
ksT, we can expand the hyperbolic function to find (b z ) oc M, a result confirmed over wide 
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range of temperatures by a MC simulation of the double-exchange model by Calderon and 
Brey. [6]. This M-dependent gauge field produces the AHE. Combining the data in Fig. 
2 with M — H curves with if//[100], we separated Rs and Ro by fitting to Eq. (1). The 
temperature dependent Rs is shown in Fig. 3. Qualitatively, this result shows very good 
agreement with Ye et al.'s prediction as well as Calderon and Brey's calculation; that is, Rs 
grows exponentially at low temperature, has an inflection point near Tc, and then shows 
a maximum at ~ Tc + 30 K before decreasing gradually at high temperature. Expanding 
Eq.(2) in the low temperature regime, we can simply describe Rs by 

Rs - -% = --T^^M^o cx exp (—E c /k B T)/T. (3) 
fx M it k B Tg/j, B /j, 

We fit the data shown in Fig. 3 by using Eq. (3) over a temperature region where the 
independent string picture should be valid, obtaining an exponential increase with Ec ~ 
1100 K ~ 3To Monte Carlo calculations for DE or Heisenberg models give different values 
for E c , but agree on the topological defect pair density at T c , (n)\ T =T c ~ 0-05 [3, 4, 6], 
fixing the exponential prefactor. With n e i « 0.5, we obtain a large, but still plausible, value 
for the spin-orbit coupling A so ~ 30 K. Our Eq is substantially smaller than expected from 
combining Ostlund's [19] analytic value An J for the Heisenberg model with the MC critical 
temperature Tc = 1.5J to get Ec = 8.4To; it is also smaller than obtained by the MC 
simulation of a classical double- exchange ferromagnet (Ec ~ 7Tc). [6] However, the energy 
for small deviations in the angle between neighboring spins is proportional 1 — cos 9 « 9 2 /2 
in the Heisenberg case, while the double- exchange energy varies as 1 — cos(#/2) ~ # 2 /8, 
decreasing the core energy to nJ for a double exchange system [20]. Using the transition 
temperature Tc = 1.06J from Calderon and Brey, we see that Ec = 3To is the expected 
core energy. 

The MC simulations show that the exponential rise in (n) gives way to an inflection point 
near Tc, causing d(n)/dT, and therefore dRs/dT, to exhibit a peak. Lau and Dasgupta 
indicated that the peak behaves ~ \t\~^, with ip w 0.65, while HJ's simulation suggested 
the weaker singularity ~ \t\~ a , where a is the specific heat exponent. Ye, et al. endorsed 
HJ's suggestion, but expected this weak anomaly to be obscured by the finite field required 
to observe the AHE[10]. We do see an inflection-like anomaly at Tc but cannot extract the 
critical behavior of the topological defects from Fig. 3. Rather, we extend HJ's scaling ansatz 
by including the scaled magnetic field and seek to collapse the Hall data as commonly done 



for other properties near a critical point. Scaling in favor of the magnetic field by setting 
Lfis/u h = 1 we obtain 

(n) = (n) reg - L^-^' v f{L x ' v t, L^'h) 

= (n) re9 - h^- a)/ps f(t/h 1/ps ,l), (4) 

where t and h are reduced T and H, respectively, f(x,y) is a scaling function of the free 
energy [21] and (n) reg is the non-critical part of the defect-pair density. Because p^ y is 
proportional to the magnetization, we eliminate the explicit field dependence in Eq.(4) via 
the usual scaling behavior of the magnetization, m(t,h) = h}/ 5 g{t/h}^ 5 , 1) to obtain the 
following scaling equation for (n) , 

(n) = (nr 9 -m^i^, (5) 

where the scaled temperature is x = t/h 1 ^ 5 . Combining this with Eq.(3), we can express 
the anomalous Hall resistivity as a function of m and scaled temperature, 

P& T = p° xy T c m(l - D{x)mP- a W). (6) 

In this equation, D(x) is the fraction term in Eq. (5) divided by (n) re9 and p xy contains all 
other constants. Fig. 4 shows the various isotherms of p xy T as a function of the reduced 
magnetization, m = M/M . At low temperature, both m and p xy increase rapidly as 
the domains are swept out. The slight positive increase in the Hall resistivity close to 
saturation is from the OHE. As temperature increases, the anomalous term dominates and 
the signature of the ordinary term becomes difficult to detect. Close to the Tq = 384 K, 
the p xy T curves reach a maximum at m ~ 0.6. The extremum of Eq. (6) is at m = 
(P/(D(1 -a + Z?)))^ 1 -^ w (0.24/D) - 316 for a = -0.17 and /3 = 0.37. The curve of Fig. 
4 follows the experimental data near the critical temperature for D(x = 0) = 1.2. The 
dashed curve represents the best fit when the LD exponent i[) = 0.65 is used in place of a; 
D = 0.8 in that case. Along x = 0, we expect p^ y — > as m — > 1 so that D(0) w 1 is not 
surprising. It is remarkable that the wide range of fields and temperatures that comprise 
the plots in Fig. 4 collapse as they do. This result strongly supports HJ's suggestion that 
d{n)/dT has the same exponent as the heat capacity and adds credance to LD's assertion 
that the topological spin defects play crucial role in the phase transition, . 
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We conclude that the AHE of CrC-2, a self-doped DE system, is consistent with Ye et 
al.'s theory and Calderon et al.'s MC simulation in which topological spin defects produce 
a strong gauge field, proportional to the magnetization and added to the applied physical 
field. Further an extended scaling hypothesis for the topological spin-defect density explains 
the behavior of p along the critical isotherm and beyond. We have argued that the AHE 
is the first experimental method able to probe the existence and properties of topological 
spin defects in a 3D Heisenb erg-like ferromagnet. The pronounced Berry phase induced by 
spin texture in a double exchange magnet has opened the opportunity, for the first time, to 
explore the fundamental role topological defects play in the critical behavior of 3D systems. 
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FIG. 1: Resistivity of C1O2 for 6— and c— directions as a function of temperature. The Psook/Piok 
ratio is 20 for b— axis and 66 for c— axis. The upper inset shows the temperature derivative of p xx 
for both directions. The rounded region is excluded as not in the critical region for fitting. The 
lower inset is a Hall pattern for this experiment to measure both 6— and c— directions at once. 



9 



I// b-axis 

o I// c-axis 




FIG. 2: Field dependence of Hall resistivity p xy (H) at various temperatures. At low temperature, 
the anomalous part cannot be observed but a slight kink is seen at low field, reflecting the compli- 
cated multiband structure. Around Tq, the anomalous Hall component is dominant and has a sign 
opposite to the ordinary Hall coefficient. The solid lines and open circles are the Hall resistivity 
for b— and c— directions, respectively. They are identical within the experimental error. 
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FIG. 3: Anomalous Hall coefficient Rs obtained by fitting of Eq. (1). Because of the difficulty in 
separating Rs from Rq, the data points in the paramagnetic regime are uncertain, the points below 
400 K are easily decomposed. The solid line is a fitted result of Eq. (3). Rs = Aexp(—Ec/kBT)/T 
is a practical fitting function. The activation energy of the topological defects and a prefactor of 
the fitting function are obtained as E c = 1103 ± 28 K and A = -2.90 ± 0.25 x 10~ 4 m 3 K/C, 
respectively. Inset: Cartoon of a pair of topological spin defects (skyrmion string) with Q = ±1. 
The hedgehog and anti- hedgehog are strongly bound, one lattice spacing apart. 
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FIG. 4: p xy T vs reduced magnetization m at various temperatures. In the critical regime, all points 
collapse onto a single curve which has an extremum atmw 0.6. The thick solid line is the best fit 
for the modified scaling function of the anomalous Hall resistivity, Eq. (6) at the critical isotherm 
(x = 0), and the obtained parameters, p® and -D(O) are -3.57 p£l cm, and 1.2, respectively. The 
dashed curve is the best fit using ip = 0.65 in place of a with D(0) = 0.8 and pIL = —5.34 fiQ cm 
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